Revivals, collapses and magnetic-pulse generation in quantum rings 
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Using a microscopic theory based on the density matrix formaUsm we investigate quantum revivals 
and collapses of the charge polarization and charge current dynamics in mesoscopic rings driven by 
short asymmetric electromagnetic pulses. The collapsed state is utilized for sub-picosecond switching 
of the current and associated magnetization, enabling thus the generation of pulsed magnetic fields 
with a tunable time structure and shape asymmetry which provides a new tool to study ultrafast 
spin-dynamics and ratchet-based effects. 

PACS numbers: 73.23.-b,78.67.-n,72.15.Lh,42.65.Re 



Introduction.- The spectacular advances in the design 
and tunability of the time-structure, amplitude, phase, 
and the shape of electromagnetic pulses 0, IS S IE 
have opened new avenues for their utilization in funda- 
mental and applied research. Inducing and monitoring 
non-equilibrium states as they build up and decay 0, Q 
are just but one example of recent applications. This 
progress in laser physics is paralleled with equally im- 
pressive development in the fabrication and manipula- 
tion of meso and nano-scale electronic systems in general, 
and of those with rin g-co nfining geometry in particular 
ES m El ES Ellli- Subjecting quantum rings to 
strongly asymmetric pulses results in the formation of 
time-dependent charge polarization and charge currents 
[iH Il7j . How these non-equilibrium phases collapse, re- 
vive and eventually relax due to coupling to other degrees 
of freedom have not yet been addressed. We demonstrate 
that, once the evolution dynamics of the pulse- induced 
charge polarization is known, a scheme can be developed 
which allows to abruptly switch on and off the charge 
current in the ring, enabling thus a fine tuning of the 
time structure and the shape of the light-induced current. 
This is important in so far as the current in the ring is 
associated with a magnetic field, i.e. a controllability of 
the current renders the generation of short (down to pi- 
coseconds) magnetic pulses. The pulse-induced magnetic 
field can be utilized to study and manipulate locally and 
in a non-invasive way the picosecond spin-dynamics [l^ 
of magnetic samples ^E HE 1^ US- In addition, we 
demonstrate how shape asymmetric magnetic pulses can 
be generated and controlled revealing thus the potential 
of our pulses for inducing a ratchet effect without spa- 
tially asymmetric potential ji^l- The simplicity of the 
proposed setup and the local tunability of the magnetic 
pulse shape and duration make our scheme a valuable 
addition to currently known accelerator-based methods 
EE I2E 1^ and those based on the generation of cur- 
rents in two photo switches excited by femtosecond laser 
pulses 0. 

For inducing non-equilibrium states so-called electro- 
magnetic half-cycle pulses (HCPs) |EIEQ are employed. 
An HCP is a highly asymmetric monocycle pulse whose 



electric field consists of a strong and short (duration Td > 
ps) half cycle followed by a much weaker but longer op- 
posite polarity half cycle. Collapses and revivals of the 
HCPs-induced non-equilibrium charge polarization, the 
charge current generation and suppression, and relax- 
ation effects are studied using the density-matrix for- 
malism IE EE ■ are interested in non-destructive 
processes, weak HCPs (few V/cm) arc applied. Thus, 
for low temperatures, T , the evolution dynamics is gov- 
erned by coupling to longitudinal-acoustic (LA) phonons. 
Electron-electron interaction effects arc suppresse d by 
Pauli blocking and energy-conservation restrictions [2y|. 
On the other hand, as we are interested in weak driv- 
ing fields optical phonons do not influence the evolution 
dynamics because their energy is well above our highest 
excited electron energy. 

Theory.- We consider TV electrons in an isolated ring 
with a mean radius rg, height zqj and width d (zq <C d <C 
To) at low T. We specifically address the case where d and 
zq are significantly smaller than the Fermi wave length of 
the charge carriers. Hence, only the lowest ground-state 
radial subbands are populated, a situation which is ex- 
perimentally feasible for semiconductor-based rings |3ll |. 
The single-particle energies associated with the angular 
motion (characterized by the quantum number m) read 
£m = fi^TO^/(2m*rQ), where m* is the effective mass. The 
particles' angular dynamics is governed by the Hamilto- 
nian iJtot = 
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Here 



is the free carrier Hamiltonian and a]^ (flm) 



phon 



E 

denote the creation (annihilation) operators. iJ, 
Yliq^^qi^q^q^^ the frce-phouou Hamiltonian, uj^ 

is the frequency of a phonon with a momentum (f, and 6t 
(6g-) are phonon creation (annihilation) operators. The 
carrier-phonon coupling is dictated by the Hamiltonian 
Hp = E,~,™,™'5rVL«m-™' +h.c., where g^-^' = 
(yr^""' Jd"^r ijit, {r )e^^^ '4^m' jr) is the electron-phonon cou- 
pling constant [25| and tpm{r) are wave functions of 
the ring carriers in the lowest radial subband. The 
bulk electron-LA-phonon coupling constant 5^""^ is de- 
termined by Iff^"""!^ = h\D\'^q/{2Vc^j,p,), where D is the 
deformation constant, cla is the LA- velocity of sound. 



2 



and ps is the lattice density. The interaction with a time- 
dependent linearly polarized external electric field Fit) 
is given hy V = -eF(t)ro X)™,™- (™l cos 0|m')oj„a™' . (j) 
identifies the electron angular position with respect to 
the pulse polarization axis and e is the electron charge. 

A key aspect of this work is that the system is ex- 
cited with a pulse which is very short on the relevant 
relaxation time scale. Hence, the excitation process can 
be separated from relaxation and dephasing which occur 
then in a field-free manner and can be monitored by mea- 
suring the emission spectrum. The excitation dynamics 
of the system is governed by Heisenberg's equation of 
motion (EOM) for the density operator Pm,m' = O'ti'-'-m', 



i.e., ihdtPrnu 



jjjcarr _|_ V , Pmi,m2\- In the interaction 



representation (p^* = Pmurn-ie 



m)*) we find 



ifidtp^l 



leF{t)ro{p^l,^^+,ci^^- 



-e„2 + i)t 



-1-0™* iCR 

f^mi + l,m2 
Pmi— l,m2 



(1) 



If the applied pulse duration Td is much shorter than 
h/{em+i ~ £m); where Sm is in the vicinity of the Fermi 
level (Ep) within a range determined by T, then the ex- 
ponential factors in Eq. hardly vary on the time scale 
Td- For low T this situation is equivalent to <C 
where r^^ ~ Sttj'o /Vp and is the carrier Fermi velocity. 
This case (called impulsive approximation (lA) [221123) 
is indeed realizable experimentally: for typical ballistic 
rings can be tens of picoseconds, while pulses down to 
subpicosecond time scale can be generated with contem- 
porary techniques 01 • Under lA, Eq. reduces in the 
basis-free form to ihdtp^^^ = [— ei^(i)ro cos 0, p™*]. If the 
duration of the negative HCP tail is much longer than 
Tp the action of this tail on the system is averaged to 
zero. Applying the pulse at t = we operate again in 
the Schrodinger picture and derive for the density matrix 
the relation (for t > Td) 



mm' 
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'm',m7 



(2) 



where p'^^i „i is the density matrix before the appli- 



cation of the pulse. C„ 



(me 



id cos (j) 



mi) 



J,„j_m(a), Ji{x) are Bessel functions, and a — 
rop/h, where p — — e J^'' F(t)dt is the momentum trans- 
ferred to the charge carriers. If for t < the system 
is in an equilibrium state with an initial density ma- 
trix P°„/,„, = /^(?7, T)6,n;m, where /^(t?, T) is the Fermi- 
Dirac distribution function [r] and T denote the chemical 
potential and the temperature, respectively) • For 
our isolated rings the chemical potential at a given tem- 
perature is uniquely determined by the number of par- 
ticles in the ring N . The pulse electric field is chosen 



as F{t) = i^o: 



exp I 



The pa- 



rameters Fq, To, and h determine respectively the ampli- 
tude, the duration, and the asymmetry of the pulse. The 
duration Td of the positive half cycle is determined by 
Td = To(l + Vl + 262 lnb'^)/b. For & = 8 the pulse has 
basically the experimentally observed shape and ratio of 
13:1 between the maximum field values of the positive 
and negative polarity parts The peak field strength 
is Fp « 0.593Fo. When rings are irradiated with these 
pulses, a charge polarization builds up which is character- 
ized by the dipole moment fl{t) = TT[erp{t)]. The com- 
ponents of the dipole moment along and perpendicular 
to the pulse polarization arc (t) ~ ero J2m ^'^[Pm+i,rn] 
and p±{t) = ero X]mlm[/0m+i,m], respectively. A de- 
tailed comparison between the dynamics of the dipole 
moment calculated using the lA and the exact numerical 
solution for different Td but fixed a endorsed our expec- 
tation that lA is well justified if Td is smaller than a 
quarter period of the dipole moment oscillations (a pe- 
riod of oscillation is determined by the energy difference 
between levels near Ep) and the HCP tail duration is 
longer than a quarter period of oscillations. Thanks to 
the lA a single parameter, the kick strength a, is suffi- 
cient for characterizing the coupling of the HCP to the 
electronic system. 

The non-equilibrium dynamics of the induced dipole 
moment is inferred from EOM for Pmi,m2 (including scat- 
tering from phonons), i.e., |25l | 



ihdtPr, 



= (er, 



E 



^l)P^ 



-mi m.3,m2 



<9T 



Using the Markov approximation and neglecting po- 
laron corrections the phonon-assisted density matrices 
^mi,m2 ^ (6gaJ,j^a,„2) - (bg) (al^Orn^) in absence of an 
external mechanical action read 12 



mi ,7712 

^9 



E 



^ [(^(f f )Pmi ,777.4Pm3 ,m2 ^gPms ,m2 Pmi ,7774] : 

where Pm,m' = Sjn,m' - Pm,m' ■ For a < 1, i.e. for weak 
excitations, we write Pmi,m2 ~ fmt5rm,m2 + Pmi,m2- As- 
suming thermal equilibrium for phonons [i.e., = nP{q), 
where n'^{q) is the Bosc-Einstcin distribution function], 
performing the sums over q, and assuming the Debye 
model for the LA-phonon spectrum, we arrive at the 
equation for our numerical calculations of the relaxation 
dynamics 



dtP7, 
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FIG. 1: (a) Dynamics of the dipole moment fi for different 
temperatures. Kick strength a = 0.2. (b) Emission spectra. 
The curves are offset vertically for clarity, (c) Dependence 
of the relaxation time of the dipole moment on temperature 
T and the number A*' of electrons, (d) Numerically calcu- 
lated current decay (dotted line) after excitation by the jTci- 
delayed sequence of two mutually perpendicular HCPs and 
current suppression (full line) by application of the second 
■jTci-delayed sequence . Duration of each HCP is = 3 ps; 
T=l K. The arrow in (a) marks the time moment when the 
stopping sequence of HCPs is applied. Parameters of the ring; 
ro = 1.35 fim, d = 50 nm, and N = 400. 



where = hcl^ps(frQ/\D\'^, 

RZ^TT:{<iZd) [n°(e)+i-r„], ee(-^,o); (4) 



otherwise; 



' ^m')/ (fic-^j^) and uj^ is the Debye frequency; 



1^/2) 



■7 if sgn(m) 



^miy) '— y JO "^-^ ^i_^2/y2 i:2[2:2-(27r)2]2 

sgn(m') and J^™ (y) ~ otherwise. 

Numerical results.- The short-pulse induced dynamics 
occurs within the following scenario. Starting from the 
ground state a short pulse drives the system to the ex- 
cited state determined by Eq. (01. This non-equilibrium 
state serves as the initial condition for the numerical 
propagation of equation (jS)), which in turn govern the 
post-pulse charge dynamics. 

For the calculations the following ring material param- 
eters (corresponding to n-GaAs) are employed; cla = 
4.79xl05cm/s, = 5.32 g/cm"^ \D\ = 8.6 eV, hu^ = 
30 meV, m* = 0.067mo (too is the free electron mass). 
Fig. n (a) shows the nanosecond-dynamics of the dipole 
moment /.j for N = 400 and different T (the oscillatory 
behavior on the much shorter time scale Tp is not resolved 
[l6|). For low T (= 0.1 K) only states with two possible 



energy values near the Ep are populated. The long-time 
behavior exhibits thus beatings and hence two peaks in 
the emission spectrum (Fig. ^(b)) are observable. With 
increasing T more levels near Ep are populated and the 
dynamics of the dipole moment shows alternating col- 
lapsed states and quantum revivals. The revival time is 
given by T^ev = 47r/i/|9^e„i/9m^| with e„i in the proxim- 
ity of Ep. The fast (classical) oscillations have a period 
Tci = 2Trh/\dem/dm\ « Ty|23- The decay of the revival 
peak values is due to the relaxation. For an estimate of 
the dipole-moment relaxation time, Trei, we utilize the de- 
cay dynamics and use the two first consecutive envelope 
maxima to determine Tj-ci from a fit by ^exp(— t/rrei). 
Fig. ^c) displays the dependence of Trei on T and on the 
number of electrons TV in the ring. Tj-ei increases rapidly 
with decreasing T and, as expected, Trei tends to infin- 
ity for T — > 0. Fig. n] (b) evidences that by measuring 
the emission spectrum the temperature-dependent decay 
dynamics can be traced. The different decay dynamics 
of the dipole moment and of the current (which can be 
monitored separately) allows insights into different parts 
of the density matrix. The dipole moment (current) is 
namely determined by the diagonal (near-diagonal) ele- 
ments of the density matrix. 

Magnetic field tuning.- Applying two time-delayed mu- 
tually perpendicular HCPs generates charge current / 
and an associated magnetization M in the ring • The 
utilization of the pulse-induced magnetic field to moni- 
tor locally, in a pump-probe manner, the ultra-fast spin 
dynamics in nanostructures or to induce a fast ratchet 
effect require a fine control of the fast switching behav- 
ior of the pulse-triggcrcd magnetic field, a task which is 
tackled here for the first time by utilizing the collapsed 
states of the dipole moments. Two weak perpendicular 
pulses with strengths ay and aj_ and delay time r (cho- 
sen as T = 21:1/4) initiate along the respective field axis 
two independent polarization dynamics /i||(i) and /ij^(i). 
Using Eq. Q we derive that the charge current change 
AI by application of the second HCP with kick strength 
a± at t = tsw is given by 



AI = -a± 



Mil ft = ^sw) 



ero 



/o, Io^eh/{m*rl). (5) 



Correspondingly, the ring magnetization is shifted by 
AM = TrrgA/, and the current-generated magnetic flux 
by A$ = jiiQroAI (fiQ is the magnetic permeability 
of surrounding material). Eq. jSJ can be exploited to 
turn the current on and off producing thus short mag- 
netic pulses, as follows; Having produced a current with 
a HCPs pair we subject (after the time Tp_p) the col- 
lapsed state to a second sequence of two pulses. The first 
pulse of the second sequence has the strength a' = q;|| 
and triggers the same dynamics of the dipole moment 
as does the first HCP of the first HCPs pair because 
the collapsed state is unpolarized. The second pulse of 
the second sequence has strength a', and is applied with 
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Numerical calculations with realistic pulse and material 
parameters are performed endorsing the feasibility of the 
predicted effects with nowadays technology. 
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FIG. 2: Magnetization pulses generated in the ring by appli- 
cation of two HOP pairs delayed by rp_p — 40 ps in (a), 
Tp-p — 120 ps in (b), and rp_p = 200 ps in (c). Kick 
strengths are a|| — a± — a]| — 0.4, a'^ = —0.39 in (a), 
a'j_ = —0.35 in (b), and a'j_ = —0.32 in (c). Magnetization 
M is normalized to Mo — -rtqIo. Dotted lines show the re- 
sult without applying the second HCP pair. Parameters of 
the calculation; ro = 0.3 /im, d = 20 nm, N=160, T=20 K, 
Td = 0.5 ps. Figures (d) and (e) show periodically alternat- 
ing time-asymmetric magnetization generated in the ring by 
application of series of sequences of HCPs with a time delay 
between sequences 1.4 ns. In (d) each sequence contains four 
HCPs having kick strengths a\\ — 0.2, = 0.2, qJ| = 0.2, 
and aY — —0.1, respectively. In (e) each sequence contains 
two HCPs having kick strengths Qy = 0.2 and a± — 0.2, re- 
spectively, is the offset magnetization associated with 
an external magnetic field Wi . Parameters of the ring and 
duration of HCPs as in Fig. [T] T = 4 K. 



polarity opposite to the second HCP of the first pair. 
|q;'j_| is chosen smaller than aj_ by a factor determined 
by the current decay during Tp_p and in such a way that 
the currents triggered by the two HCP pairs cancel each 
other. Fig. ^c) demonstrates numerically the accuracy 
of this scheme of current switching. The abrupt current 
switching opens the way for the tunability of the time- 
structure of the magnetic pulses. The magnetic pulses 
duration is stretchable depending on Tp_p (cf. Fig. EJa)- 
(c)): As inferred from Figs.|Il;d)[ll;a)-(c), if Tp_p < t^^'', 
where t^^^'^ is the time constant for the current initial 



decay dynamics, almost rectangular ps magnetic pulses 
are produced. For rp_p > r^™'' asymmetric pulses are 
triggered, which is of importance when studying ratchet 
effects |2J|. Fig. |21 demonstrates the level of tunability of 
the pulse shape and duration by varying experimentally 
accessible parameters. The magnitude of the triggered 
pulses is enlarged for smaller rings and/or stronger kicks. 
Also an appropriate arrangement of a collection of rings 
allows additional tuning of the magnetic pulses. 

In summary, we demonstrated how relaxation, col- 
lapses, and revivals in quantum rings can be studied by 
means of electromagnetic pulses and how these phenom- 
ena can be exploited to generate with current technology 
magnetic pulses with tunable time and shape structure. 
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Many-body effects on the ground-state distribution func- 
tion are suppressed due to Pauli's principle. 



